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I. INTRODUCTION 

The existence of massive neutrinos provides concrete evidence for physics beyond the Standard Model (SM). 
Similarly, the explanation of observed galactic rotation curves in terms of gravitating dark matter (DM) further 
suggests the SM is incomplete. Efforts to explain these two key evidences for new physics are varied, though 
an interesting approach is to seek a common or unified framework that simultaneously solves both puzzles. 
For example, if small neutrino masses are realized via radiative effects [I|, it is conceivable that DM plays a 
role in generating the masses, allowing a type of unihed description for massive neutrinos and DM. This is the 
motivation for the models of Krauss, Nasri and Trodden (KNT) and Ma Q . Both models extend the 
SM so that neutrino masses are generated radiatively with DM propagating in the loop diagram. In order to 
ensure DM stability (and preclude tree-level neutrino mass) a Z 2 symmetry is also imposed. 

There are a number of generalizations of this basic idea which similarly extend the SM to allow radiative 
neutrino mass via couplings to DM 1M3- In common with the KNT and Ma models, the generalized models 
also require the imposition of a new symmetry to render the DM stable.^ However, it is interesting to consider 
models where DM stability instead results from an accidental symmetry, in accordance with our experience 
from the SM, where proton stability manifests the accidental baryon number symmetry. 

In this work we present a model of radiative neutrino mass that automatically contains an accidental Z 2 
symmetry and thus admits a stable DM candidate. The model realizes a simple unified framework for the origin 
of neutrino mass and DM while imposing only a minimal symmetry structure, namely that of the SM. Neutrino 
mass appears at the three-loop level via a diagram with the same topology as the KNT model, while the DM is 
a neutral fermion with a non-trivial charge under the accidental Z 2 symmetry. The model requires heavy DM 
(Mdm ^ 20 TeV) and may be probed via DM direct-detection experiments and future /r —>■ e -I- 7 searches. It 
also predicts a charged scalar that can appear at the TeV scale. 

The layout of this paper is as follows. The model is introduced in Section|TTl We calculate neutrino masses and 
discuss important constraints in Section IIIIl Relevant information regarding the DM is discussed in Section IIVI 
while our main numerical analysis and results appear in Section [V] Conclusions are drawn in Section IVTl 
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^ In some models the DM is merely sufficiently long-lived, rather than absolutely stable. This does not require a new symmetr y bu t 
instead relies on technically-natural parameter hierarchies (either among mass parameters [T^ or dimensionless couplings flOH '). 
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II. THE MODEL 
A. Field Content 


We extend the SM to include a charged singlet scalar, 5+ ^ (1,1,2), a scalar septuple!, cj) ~ (1, 7, 2), and three 
real septuple! fermions, JFi ~ (1,7,0), where i = 1,2,3, labels generations. We adopt the symmetric-matrix 
notation for the septuplets, writing the scalar as (j^abcdef, with a,b,... € {1, 2}. The components are given by 
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where and (j) are distinct fields, (j) ^ , and similarly (f> ^ For the septuple! fermions, 

denoted as Tabcdef, we have 
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The superscript “c” denotes charge conjugation and the numerical factors ensure the kinetic terms are canoni¬ 
cally normalized. With these helds, the Lagrangian contains the terms 


3 -^SM ~ 2 7?7? {dia 't’^aR faf) Lfj S'^ + H.C.} — V{H,S,(I)), 


( 3 ) 


where lepton flavors are labeled by lower-case Greek letters, a, f3 £ {e, r}, and L (en.) is a SM lepton doublet 

(singlet). The scalar potential is denoted as V{H,S,(j)). Note that the exotics (j) and T do not couple directly 
to the SM neutrinos, though they shall play a key role in generating neutrino mass. 

The explicit expansion for the fermion mass term is: 


{y!)abcdef My {J^j)p<irstu e"’ e'”' 7“ + H.c. 

= -M,,|7f^F-++++:^F-+++:^F-+ + iFfF-o|, (4) 

where we defined: 

.F+++=F-+++ + F-+++, F-++=F-++-F-++, F-+=F-++F-+, j-o = jro _ (5) 

Clearly is a Majorana fermion, while the other six components of F" partner-up to give three massive charged 
fermions (per generation). Without loss of generality, we choose a diagonal basis for the fermions, such that 
Mij = diag(Mi, M 2 , M 3 ), with the masses ordered as Mi < M 2 < M 3 . We shall see below that F does not 
mix with the SM leptons, to all orders of perturbation theory, so Eq. (0 describes the mass eigenstates, which 
should be used in the Yukawa terms in Eq. ©. The lightest neutral fermion will play the role of DM, and we 
denote its mass as = Mi. 


B. An Accidental Symmetry 


The model contains an exact accidental Z 2 symmetry with action: 


{cjb, F} ^ {-<)., -F}. 


To see this, note that the potential can be written as 


( 6 ) 


V{H, S, cjb) = V{H) + V{cjb) + V{S) + Vm{H, S) + Vm{H, </.) + Vm{S, </.). 


( 7 ) 
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The first four terms trivially preserve the discrete symmetry, while the explicit forms for the last two mixing 
potentials are^ 


Vm{H, Cj)) = XhMHT' + XhMHT’ ^a'bcdef, (») 


and 


Vm{S, = Xs^SWrr^^^'^^^abcdef + ^ (5" 0abcde/</'a'b'c'd'e'/'e“'e“'e“'6^^' + H.C. (9) 

These potentials also preserve the symmetry defined by Eq. dfi]) . Note that there appears to be a third distinct 
way to contract the SU{2) indices in the mixing potential Vra{S, tj))^ namely 


o—( j.^XQ'bcdef 1 1 a'd' h'e'c' f 

^ (0 / Yabca'b'c'4^defd'e'f'^ ^ ^ 


( 10 ) 


This would explicitly break the Z 2 symmetry. However, this term is odd under the simultaneous interchange of 
the sets of dummy indices {a,b,c} O {d, e,/} and {a',b',c'} O {d',e',/'} [l^, and thus vanishes identically. 
The full theory therefore preserves the accidental Z 2 symmetry defined by Eq. ([HI) and the model automatically 
contains an absolutely stable particle that is a DM candidate. The Z 2 symmetry also prevents mixing between 
T and the SM leptons. To the best of our knowledge this is the first such model of radiative neutrino mass with 
an accidental symmetry that automatically gives a DM candidate. 

At tree-level the components of J- are mass-degenerate, while the components of 4> experience a mild splitting 
due to the A// 02 -term in Vm{H, cj)). Eor > O(TeV) this mass-splitting is not significant and is essentially 
negligible for Xh<i ,2 ^0-1. Thus, to good approximation the components of Ti are degenerate at tree-level, with 
masses Mi, as are the components of (j) (with masses M^). Radiative corrections remove these mass degeneracies; 
loops containing SM gauge bosons give small mass-splittings for the components of T, leaving as the lightest 
exotic fermion. Similar splittings are induced for the components of 4> which are readily calculated with the 
results of Ref. [l^. For most purposes in this work these tiny splittings can be ignored. 

The model contains two distinct Z 2 -odd DM candidates, namely and cffi. However, cffi has degenerate real 
and imaginary components and also couples to the Z boson. This leads to tree-level Z boson exchanges that 
are incompatible with direct detection constraints. Thus, can be excluded as a DM candidate, leaving 
as the sole DM candidate in the model and restricting one to the parameter space with = Mi < M^. The 
SM Higgs develops a nonzero vacuum value, (H) 7 ^ 0, breaking the electroweak symmetry in the usual way. 
Furthermore, in the parameter space with (cj)) = 0, which preserves the discrete symmetry, the p-parameter 
retains its standard tree-level value.^ 

Note that the inclusion of non-renormalizable dimension 5 operators [D = 5) can break the accidental 
symmetry. In particular, the D = 5 operator HH(j>^(j>^(j> would allow the DM to decay. This feature is not 
specific to the present model; one expects global symmetries to be broken by gravitational effects [l^ so non- 
renormalizable operators will, in general, break global symmetries. This is true even in related models which 
impose e.g. a global Z 2 or 17(1) symmetry. In our model, the fate of the accidental Z 2 symmetry is analogous 
to the fate of the accidental baryon number symmetry in the SM. The latter is broken by D = 6 operators, 
leading to proton decay. However, proton longevity can be ensured by the details of the UV completion, giving 
either a long-lived or absolutely stable proton. Unsurprisingly, the situation is similar for our DM candidate. 
In this work we focus on the calculable renormalizable theory. 

We note that a number of earlier works have studied larger multiplets in connection with neutrino mass [T^ . 
MM (for related phenomenology see Ref. M)- In particular, Ref. [l^ recently considered stable quintuplet 
fermionic DM in a three-loop model of neutrino mass.^ 


III. 

The 


THREE-LOOP NEUTRINO MASS AND LEPTON FLAVOR VIOLATING CONSTRAINTS 

combination of the Yukawa Lagrangian and the terms 

V{H,S,4>) D +H.C. 


= y(5-)^{^++++r- 


A-l— 


■k++jfi 


b+A+ 


} + H.c. 


( 11 ) 


^ The second term is equivalent to the standard {H^riH){(l)^Ti(p) term, where and Tj denote SU{2) generators for the distinct 
representations. 

^ We shall see below that the septuplets must be heavier than the TeV scale; given the very small mass-splittings, relative to the 
weak scale, this should ensure that the new contributions to the oblique parameters are negligible. Also, similar to other models 
with large multiplets, the SU{2)]^ coupling constant encounters a Landau pole in the UV, due to the heavy septuplets. 

^ Interestingly, the model of Ref. [T3l gives an accidental Z 2 symmetry after imposing a separate Z 2 symmetry. 
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FIG. 1: Diagram for radiative neutrino mass, where S and (f> are new scalars and T is an exotic fermion. The lightest 
component of is a stable dark matter candidate. 


in the scalar potential, are sufficient to explicitly break lepton number symmetry. Consequently SM neutrinos 
are Majorana particles that acquire radiative masses at the three-loop level, as shown in Figure [1] In the limit 
where the mass-splittings among components of (j) and are neglected, the calculation of the loop-diagram 
gives 




7As rrijins 
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( 12 ) 


where the function F encodes the loop integrals Q and Mg is the charged-singlet mass. 

Neutrino masses calculated via Eq. m must satisfy the data from neutrino oscillation experiments and 
reproduce the following best-fit regions for the mixing angles and mass-squared differences: 5^2 = 0-320^g Qi[ 7 , 
sia = 0.43t°:°|, = 0.025lHo!> Ami, = 7.62l0;}9 x 10-5eV^ and \Aml^\ = 2 . 551^9 x IQ-^eV^ M- 

Matching to these experimental values reveals the regions of parameter space where the model gives viable 
neutrino masses. 

The Yukawa couplings 9 ia generate flavor changing processes like ^ e -I- 7 . Calculating the corresponding 
diagrams in the limit where the mass-splittings are neglected, and including the diagram containing the singlet 
S, gives 


Bin ey) ~ 
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where F 2 {x) = [1 — 6 a; -I- 3x^ + 2x^ — 6x^ loga;]/[ 6 (a: — 1)"^] and v is the vacuum expectation value of F[. The 
related expression for B[t ^ ^ + 7 ) is obtained by a simple change of flavor labels in Eq. (fO|) . Replacing the 
final-state electrons with muons in the diagram for /a —>■ e -I- 7 gives the one-loop contributions to the muon’s 
anomalous magnetic moment. In the limit where the radiative mass-splittings are neglected these give 


5afi = 


IGtt^ 



Ml 


F2{MllMl) + Y. 



(14) 


The last term is due to the charged scalar S. 

A further constraint of (A4^)ee ^ 0.35 eV follows from null-results in searches for neutrino-less double-beta 
decay [l^, though analysis shows that this constraint is readily satisfied in the model. This constraint is 
expected to improve after next generation experiments, with an anticipated precision of {Mv)ee ^ 0.01 eV [^ . 


IV. DARK MATTER 


A. Relic Density 


As mentioned already, the only viable DM candidate in the model is the lightest neutral fermion There 
are two classes of interactions that can maintain thermal contact between the DM and the SM in the early 
universe. Interactions mediated by the scalar (j) have the cross section 


a(2J-o^£+C) 


\9iM^MIJM^^ + MI) 

4877 (M2^ + M|)4 


(15) 


where Vr is the DM relative velocity, in the centre-of-mass frame. Note that there are no s-wave annihilations 
when final-state lepton masses are neglected, as the DM is a Majorana fermion. There are no coannihilations 
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mediated by 0, though given the small radiative mass-splittings, one should include the annihilations of singly- 
charged fermions: 




\gipgla? + 

48^ 


and similarly for the higher-charged fermions 

^ £+C) = <jli = ^±' 


X Vr = 


olB 

> 


a{T—F+++^ ip-) = = af. 


(16) 


(17) 

(18) 


There are also processes mediated by SU{2)l gauge bosons, which can be calculated in the limit of an exact 
SU{2) symmetry. The corresponding cross sections can be obtained with the results of Ref. [2l|. Due to the 
small mass-splitting among the components of J-i, one should also include coannihilation processes. Adding 
annihilation and coannihilation channels together in the standard way gives (2^ 


(7eff{2J- —>■ SM) X Vr 


1 

9lff 


CTlV X Ur -t- ^ l^O O'o.o + 25 ± + 25 ±± 0-±± + 2 g±±± X Vr 


(19) 


where the mass-splittings among fermion components are neglected and the SU(2)l channels give 

In the above, geff = ffo + 2g± -K 2g±± -H 2g±±±, with go = g± = g±± = g±±± = 2. 

In principle one can calculate the mass range that gives a viable DM relic density using the above expressions. 
However, the cross section into gauge bosons may be significantly enhanced by the non-perturbative Sommerfeld 
correction (23l - l^ . One must solve the Schrodinger equation in terms of a non-relativistic bound state of two 
DM particles in order to estimate the non-perturbative Sommerfeld correction. The calculation is somewhat 
involved, though the correction has been calculated for several SU{2)l multiplets in Ref. [2l| and the effect 
is found to be important for larger multiplets. The enhancement of the cross section influences the DM mass 
required to give the observed relic density as the DM mass is the unique parameter that can control the cross 
section when the annihilation cross section is dominated by gauge interactions.® For example, the DM mass 
is shifted from 3.8 TeV to 9.5 TeV for a fermion quintuplet with Y = 0, from 5.0 TeV to 9.4 TeV for scalar 
quintuplet with Y = 0, and from 8.5 TeV to 25 TeV for scalar septuplet with V = 0 [2^. A similar enhancement 
is expected for the fermion septuplet DM with V = 0 in our model, though a detailed calculation is beyond 
the scope of this work. Guided by the results listed in Ref. [2^ we expect the Sommerfeld enhancement will 
increase the requisite DM mass by a factor of approximately 3. As we shall see, this suggests the required DM 
mass should be ^ 20 — 25 TeV when the Sommerfeld effect is taken into account. 

The DM annihilation processes which induce monochromatic gamma-rays are also enhanced by the Som¬ 
merfeld correction in the present universe. This can be a significant signature of DM as an indirect detection 
signal. Since the DM mass is predicted around =20^25 TeV in our model, after including Sommerfeld 
correction, monochromatic gamma-rays at Ej = could be detected by future gamma-ray experiments such 
as CTA 


B. Direct Detection 


There is no tree-level coupling between DM and quarks. However, W boson exchange gives three one-loop 
diagrams which can produce signals at direct-detection experiments [^. There are both spin-dependent and 
spin-independent contributions to the scattering, however, spin-dependent contributions are suppressed by the 
heavy DM mass. As we consider relatively heavy values of Mdm > 1 TeV, the spin-dependent contributions 
can be neglected. Therefore spin-independent scattering dominates and the cross section is determined by SM 
interactions: 


asi E°N) 






(21) 


^ This is expected in the present model, due to the relatively large value of crvy in Eq. iiM . 
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FIG. 2: The DM and charged scalar masses versus the scalar septuple! mass for the case with no Sommerfeld enhancement. 
The blue line at Mdm = 7.2 TeV gives the best-fit value for in the limit gta —>■ 0. 


The DM scatters from a target nucleus A of mass and the standard parametrization for the nucleon is 
adopted: 


{N\^mqqq\N) = fuiN- ( 22 ) 

Here ruN is the nucleon mass and / = J^qfg is subject to the standard QCD uncertainties. For / Ri 0.3, the 
cross section for the one-loop processes is crsi ~ 4 x 10“''’^cm^, which is just beyond the sensitivity of LUX [2^ for 
heavy DM with Mdm ~ 25 TeV. Note, however, that recent lattice simulations suggest a somewhat lower value 
of strange content fg ~ 0.043 ± 0.011 [2^, which, when combined with cancellations from two-loop diagrams, 
gives a smaller cross section of trsi Ri 4 x 10“"^®cm^ [1^.® In either case, the result is beyond the current 
sensitivity of LUX, though future discovery prospects for the DM candidate can be considered promising. 


V. NUMERICAL RESULTS AND DISCUSSION 

As already mentioned above, to determine the viable DM mass range one should include the Sommerfeld 
enhancement. However, as a first task we perform a numerical scan of the parameter space without the 
Sommerfeld enhancement, determining the favored DM mass range. We subsequently include a simple estimate 
of the effect. 

For the numerical scan we seek regions of parameter space that satisfy the previously mentioned constraints, 
while simultaneously giving neutrino masses and mixings in agreement with the experimental values and a DM 
relic density within the range flouh^ ~ 0.09 — 0.14. We consider the free parameter values 

lUpf , \g^c.\^ < 9, 500 GeV < Mdm < 10 TeV, 

100 GeV < Ms < 10 TeV, > Mdm- (23) 

The results for the values of Mdm, Mg and are shown in Figured We find that viable neutrino masses 
can be obtained for a large region of parameter space, though the DM mass should be confined to the tidy 
range of 7.18-7.31 TeV for the relic density to match the observed value. This region is somewhat tighter than 
the corresponding region for the related models with triplets @ and quintuplets 0 , due to the fact that the 
cross sections for annihilations mediated by the couplings gia, namely Eqs. (Ha-dlHl) jUre smaller compared to 
the contribution of SU{2)l gauge bosons (1^ . In the triplet and quintuplet cases @,[3 the charged lepton 
contribution is non-negligible, allowing a greater spread for the DM mass interval. 

The Sommerfeld enhancement is expected to increase the required DM mass by a factor of roughly 3. There¬ 
fore, in order to approximately take this effect into account, we redo the numerical scan with the DM mass in the 
relic density replaced by Mdm/ 3, searching for parameter space that gives viable neutrino masses and mixings 
and is consistent with low-energy constraints. This approach only provides a rough approximation for the value 
of the DM mass but, importantly, it allows us to discover if the requisite heavier values of Mdm and M^ are 
compatible with the low-energy data. Note that, because the relic density calculation has a reduced sensitivity 


We estimate the two-loop effect with a simple scaling of the results in Ref. l26l . 
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FIG. 3: Numerical results for the case where the Sommerfeld enhancement is considered. Left: the Yukawa coupling 
range that satisfies all previous constraints. The dashed line represents the case where all couplings are close in absolute 
value, i.e., mm(|/|) = max{\f\). Center: the masses Mdm and Ms versus the septuplet scalar mass M^. The blue line 
gives the Mdm best-fit value for Q,DMh? when gia 0. Right: The constraints from lepton flavor violating decays scaled 
by the experimental bounds versus the charged scalar mass. Here the muon anomalous magnetic moment is smaller than 
the experimental bound by more than one order of magnitude. 


to the couplings gia (as DM annihilations in the early universe are dominated by SU{2)l annihilations), the key 
question is whether there is viable parameter space that achieves neutrino mass and satisfies the constraints, 
given the heaviness of the DM. Our approach allows us to answer this question and a small shift in Mdm should 
not significantly affect the conclusion. 

Performing the modified numerical scan produces the new results shown in Figure [3] There is considerable 
parameter space that satisfies the constraints with the DM mass in the range 19.7-23.1 TeV, centered around 
the value of Mdm = 21.7 TeV, which is preferred in the limit gia = 0. The scalar </) must now be heavier 
than 19.9 TeV, while the charged scalar singlet S can remain as light ^ 500 GeV, similar to the case without 
Sommerfeld enhancement effect. One observes that the branching ratio B{t — >■ /r -|- 7 ) is smaller than the 
experimental bound by 4-6 orders of magnitude while the constraint of B{g, —>■ e -|- 7 ) < 5.7 x 10“^^ is more 
severe. In particular, it is evident that improved measurements of i3(/i —>■ e -I- 7 ) are capable of excluding the 
model. Though not shown in the figure, the preferred regions of parameter space are not ruled out by the data 
on the anomalous magnetic moment of the muon; the extra contribution from the exotics can contribute to the 
observed discrepancy, though it cannot explain it entirely [HI, |32| . 

We note that with only two generations of fermions JY (gsa = 0), the bound on —>■ e -|- 7 ) is violated. 

Therefore three generations of Jy are required to remain consistent with constraints from lepton flavor violating 
processes. Also, the neutrino data prefers that one does not introduce large hierarchies between Mdm and the 
other exotic masses, M 2,3 and M^, with M^^ 2,3 ~ 0(1 —10) x Mdm preferred. The exotics are therefore clustered 
near Mdm- Finally, we emphasize that the preferred values of Mdm should only be taken as a guide, though our 
analysis clearly shows that one can satisfy the low-energy constraints with the required heavier values of Mdm 
and Mff,. 


VI. CONCLUSION 

We presented an original model of radiative neutrino mass that automatically contains an accidental Z 2 
symmetry and thus provides a stable DM candidate. This gives a common description for neutrino mass and 
DM without invoking any symmetries beyond those present in the SM. The DM is the neutral component of 
a septuplet fermion ^ (1, 7, 0), and should have mass Mdm ~ 20 — 25 TeV. The model can give observable 
signals via flavor-changing leptonic decays and DM direct-detection experiments. It also predicts a charged 
scalar S that can be at the TeV scale and within reach of future colliders. 
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